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ABSTRACT 

Orbiting spacecraft such as large space antennas have to maintain a highly 
accurate shape to operate satisfactorily. Such structures require active and 
passive controls to maintain an accurate shape under a variety of disturbances. 
This paper is concerned with methods for the optimum placement of control 
actuators for correcting static deformations. In particular, attention is 
focused on the case where control locations have to be selected from a large set 
of available sites, so that integer programming methods are called for. The 
paper compares the effectiveness of three heuristic techniques for obtaining a 
near-optimal site selection. In addition the paper presents efficient 
reanalysis techniques for the rapid assessment of control effectiveness. Two 
examples are used to demonstrate the methods: a simple beam structure, and a 

55m space-truss-parabolic antenna. 







INTRODUCTION 


Z 


In the design of large space antennas, one of the most stringent design 
requirements Is that of surface accuracy [1,2], While studies have shown that 
In some cases high-surface accuracies may be maintained with passive methods 
[3il, It Is expected that for many applications active controls may be needed. 

The disturbances which affect the shape of space structures are of two types. 

One type Is transient which leaves the structure unchanged once damped out. 

Such disturbances usually call for active or passive controls which enhance the 
damping of the structure. The second type of disturbance Is typified by fixed 
deformations such as those due to manufacturing errors [4] or those which are 
slowly varying and may be considered quasi-static. These latter disturbances 
may be offset by slowly-applied, long-acting corrections. Most research to date 
has concentrated on the first type of disturbance and the use of damping 
actuators [5]. There has been less research on controlling quasi -steady 
disturbances. 

Much of the work reported on active control of quasi-steady disturbances is 
related to active control of optical systems such as mirrors (see [6] for a sur- 
vey of the state of the art of 1978). Generally, the actuators employed are 
force actuators (e.g, [7-11]). Bushnell [7] characterizes some such actuators 
(e.g. [12,13]) as displacement actuators because they are stiff enough to 
enforce a prescribed displacement at a point. Another variation of the force 
actuator In a truss structure is one which effects a change in the length of a 
member by reeling a cable in or out or by using a screw mechanism. This 
approach is used on some antennas (e.g. [14]) to correct fabrication errors, 
albeit on the ground rather than in orbit. A recently-proposed alternative [15] 
Is the use of applied temperatures on the structure. 
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The present paper describes a follow-on effort from that of [15], namely 
the optimal placement of force or temperature control actuators In a flexible 
structure to correct static surface distortion. When the sites available for 
location placement are a continuum, the problem can be treated by standard con- 
tinuous optimization techniques. For example, [16] contains a survey of tech- 
niques employed for actuator placement In vibration control, and [17] describes 
actuator placement for static shape control. In many practical problems, only 
discrete sites are available. The optimal selection of the locations becomes an 
Integer programing problem which Is usually much more difficult and costly to 
solve than a continuous optimization problem [18]. This discrete site problem 
has received relatively little attention. 

Because of the high cost of solving Integer programing problems rigorously, 
there Is merit In considering heuristic site selection techniques which obtain 
near-optimal solutions at a relatively low computation cost. In [16] and [19] 
heuristic methods were developed In connection with actuator and sensor 
placement for vibraton control problems. An Important adjunct to these 
technqiues Is the rapid evaluation of the effect of adding or eliminating 
actuators. Reference [16] describes approximate methods for this type of evalu- 
ation In vibration control problems. The present paper describes two heuristic 
algorithms for actuator placement for static shape control. The paper also 
develops rigorous yet rapid analysis methods to evaluate the effects of adding 
or deleting actuators. 

The methods discussed In the paper are applied to two problems: a free- 

free beam and a 55m space truss parabolic antenna reflector. Studies are made 
of the effect of starting points on the final design produced by the 
algorithms. Efficiency and effectiveness of the methods are evaluated and 
compared. 
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coefficient matrix in Eq, (4) 

vector added to matrix A due to added actuator (Eq. 16) 
diagonal term in matrix A corresponding to added actuator (Eq. 16) 
unit vector with unity in ith position 

ratio of corrected to original rms displacement (Eqs. (7) and (8)) 

number of available locations for actuators 

number of actuators 

right-hand-side vector in Eq. (4) 

rms value of distortion (Eq. (9)) 

vector of incremental temperatures 

vector of incremental temperatures for nominal system 

incremental temperature of ith actuator 

displacement due to a unit temperature increment in ith actuator 
rms value of displacement 
reference volume 
residual displacement (Eq. 1) 

Lagrange multiplier (Eq. 11) 
augmented function, Eq. (10) 
disturbed shape 

continuum occupied by structure 

REVIEW OF STATIC SHAPE CONTROL METHOD 

The methods discussed in this paper are applicable to both linear force and 


temperature controls. For completeness the equations for temperature control 
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from [15] are summarized herein. The reader is directed to [15] for a similar 
derivation for force controls. 

The structure is assumed to be in earth orbit and possess rigid body 
degrees of freedom. The structure is defined over some region n and it is 
assumed that its desired shape has been distorted by an amount 4*(Q) where Q 
is a point in n. It is also assumed that the distortion can be accurately 
measured. The disturbance is corrected by prescribing temperatures at n 
high-thermal-expansion inserts (actuators) placed in the structure. The 
disturbance is assumed to be slowly varying so that the actuator inputs may 
be calculated by a quasi-static analysis. 

The residual displacement is the sum of the disturbed shape and the 
correction 

s = + “i ^1 (1) 

where T-j is the change in temperature of the ith actuator with respect to 
the temperature at which 41 is measured and u-j is the displacement due to a 
unit value of T-| . 

The best values of T^ are those which most effectively nullify 41 and 
cause cu to be close to zero. A common measure of the smallness of 
is based on the rms value 



where Vq is a reference volume. The necessary condition for a minimum is 
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( 3 ) 
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Equation (3) is a system of n linear algebraic equations for the control 
temperatures and may be written as 

AT = r (4) 

where the component aij of the matrix A is 


- "i * 

-4 


(5) 


( 6 ) 


and the jth component of the right-hand-side, rj is 

Pj = -A, . ujda 

The ratio of controlled to uncontrolled rms distortion, g is given by 

, e.^dft 

9^ = - 7 ^ (7) 
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It follows from Eqs, (1), (4), (7) that 
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THE EFFECT OF REMOVING OR ADDING ONE ACTUATOR 


The optimization algorithms used in this paper always compare a given 
configuration to another which is identical except that one actuator has either 
been removed or added. To reduce the computational cost of these algorithms. 
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quick reanalysis procedures are derived below to assess performance for these 
special cases. 


Removing an Actuator 

Eq. (4) may be obtained by minimizing g^ with respect to T. Removing the 
ith actuator can be simulated by performing the minimization under the 
constraint that T-j = 0. Employing Lagrange multipliers we look for stationary 
points of 4) where 

. _ 1 ^ 2 -T - 1 .. 2„ 2 Tt /,n\ 

♦ ■ 2 '^mo^ - j r^^g - Xe^ T (10) 

where ei is a vector with unity in the ith row and zeros elsewhere. 

The conditions of stationarity of <|> are 

A T - r - Xe^ = 0 (11) 


then from Eq. (11) 


e^ T » 0 




( 12 ) 

( 13 ) 


where Tq is the nominal vector of actuator temperatures (Tq = A-lp), 
From Equations (12) and (13) it follows that 
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where aii"^ is the ith diagonal of A"^. Then from Eqs. (8), (13), and 
(14) 
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Adding an Actuator 

Adding an actuator requires Increasing the order of A. Eq, (4) becomes 



where C Is the additional column and d the new diagonal element. From the 
expansion of Eq. (16) we obtain 


T = A-l r - A-lcT„^j = T„ - A-‘ 

(17) 

T - '•n*! ■ 

d - cV^C 

(18) 

,2, 1 
^(9 ) “ T.-l 2 

(19) 

d - c'A r ^ 

mo 



SELECTION OF OPTIMAL 
ACTUATOR LOCATIONS 


The problem of selecting n actuator locations from a set of m available 
sites can be formulated and solved by standard Integer programming techniques 
[18], However, these tend to be extremely costly when the number of available 
sites Is large compared to the number of actuators. This Is due to the large 
number of possible combinations of placement configurations. For example, the 
number of possibilities for choosing 20 actuator locations from 100 sites Is 5.4 
X 10^^. This paper proposes two heuristic algorithms which Improve a trial 
set of actuator locations at a moderate computational cost. Additionally, a 
study is made of another heuristic algorithm due to DeLorenzo and Skelton ([16] 
and [19]). 



9 

The two proposed algorithms start with a configuration Iq which con- 
tains the desired number of actuators (n). The Worst-Out-Best-In (WOBI) 
algorithm first independently removes each of the n actuators to find the 
"worst" actuator i.e. the actuator which can be removed with the least detrimen- 
tal effect on performance. Then the worst actuator is moved to each of the 
locations outside Iq and tested by including it with the n-1 remaining 
actuators. The best location replaces the removed one. A total n+(m-n)=m 
configurations are analyzed in each iteration. Iterations continue until no 
improvement is possible. The Exhaustive Single Point Substitution (ESPS) algo- 
rithm is more thorough than the WOBI algorithm. In an iteration, it moves each 
actuator in turn from Iq to each of the m-n unused locations and analyzes 
performance at each trial location. The best of these n(m-n) configurations 
replaces Iq 3nd iterations continue until no improvement is possible. ESPS 
is more expensive per iteration than WOBI because n(m-n) > m. 

The DeLorenzo algorithm starts with a configuration where all m locations 
are occupied by actuators. The least effective location is found by removing 
one actuator at a time. This least effective actuator is removed and the 
process is repeated with m-1 actuators. The process is repeated m-n times 
until the number of actuators is reduced to n. The total number of 
configurations analyzed by DeLorenzo's algorithm is (m+n)(m-n+l)/2. Details of 
the three algorithms are given below. 

Worst-Out-Best-In (WOBI) Algorithm 

1. Select an initial configuration Iq of n actuators, and calculate 
the rms reduction factor gg. 

2. Calculate the rms reduction factor, g, for each of the n 
configurations of n-1 actuators obtained by removing one actuator 
from Iq* 
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3. Select the actuator which when removed has the least effect on g as 
the "worst" actuator. Remove it to produce a configuration Ij of 
n-1 actuators. 

4. Calculate g for the m-n configurations of n actuators obtained by 
placing an actuator at any of the available locations outside Iq* 

5. Label the configuration with lowest g as and the corresponding 
9 as gj. 

6. If gi >. 9o convergence Is obtained and Iq Is the best 
conf1gu7a1on. 

7. If gi < 9o ^0 “ h 9o “ 9l 9° 


Exhaustive Single Point Substitution (ESPS) Algorithm 

1. Select an Initial configuration Iq of n actuators and calculate 

9o- 

2. Calculate g for each of the (m-n)n configurations obtained by 
replacing one of the locations of !« by a location outside of 
lo- 

3. Label the configuration with lowest g as Ij and the corresponding 

9 as gj. 

4. If gi >. 9o convergence is obtained and Iq Is the best 

configuration. 

5. If gi < 9o set Iq = Ij and gQ = gj and go to step 2. 


DeLorenzo's Algorithm 

1. Calculate go an Initial configuration Iq composed of m 
actuators (an actuator at every available location). Set nj = m. 

2. Calculate g for all nx - 1 actuators obtained by removing one 
actuator from Iq. 

3. Select the best configuration and label It Iq and the corresponding 
9 as gQ. 

4. If nj - 1 = n the process is completed. 

5. If ni - 1 ^ n go to step 2. 
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Sequential Application of WOBI and ESPS 
To minimize the variation In WOBI and ESPS results due to the selection of 


the Initial configuration the following strategy Is useful. Each method Is 
first applied for a small number of actuators and then the number of actuators 
Is Increased by one until the required number Is reached. Each time, the best 
configuration obtained for n actuators Is augmented by the first available 
site and used as the initial configuration for the selection of the n+1 
locations. 

RESULTS AND DISCUSSION 
Free-Free Beam 

The first example is a free-free beam (fig. 1) Initially straight but dis- 
torted Into a shape described by a cubic polynomial. Results of using actuator 
locations obtained with the WOBI and ESPS algorithms were compared with each 
other and with the performance of a set of equispaced actuators. The number of 
actuators (n) was either six or eight. The number of available sites (m) was 
20, 40, and 80. For each combination (n,m), ten arbitrary Initial configura- 
tions (Iq) were used to assess the variations due to the initial configuration 
(see Table 1 and figure 1(b)). 

The results summarized In Table 2 reveal the following characteristics of 
the WOBI and ESPS algorithms. The WOBI algorithm Is much more sensitive to the 
Initial configuration. In fact for some initial configurations it did not do as 
well as the equispaced configuration. The scatter in the ESPS algorithm Is much 
milder, and the improvement In its performance over the equispaced solution 
ranged up to 38 percent. The number of Iterations for convergence of both algo- 
rithms was small (typically three) and as a result, the number of configurations 
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tions which had to be analyzed Is small. For example. In selecting 8 sites from 
80 available (equivalent to 4 out of 40 because of symmetry) there are 91,390 
possible combinations. Three Iterations of the UOBI technique check 120 of 
these combinations while three Iterations of the ESPS technique check 432 
combinations. The best and worst locations obtained by UOBI and ESPS are shown 
In Figure 2. The disturbed shape and typical corrected shapes are shown In 
Figure 3. 


Antenna Reflector 

The second example is a 55m space- truss parabolic antenna reflector shown 
in Figure 4. The antenna Is assumed to be constructed of graphite epo^y while 
the control elements are aluminum. The reflector is distorted from Its Ideal 
shape by thermal deformation due to orbital heating. The temperature history of 
the lower and upper surfaes of the antenna is shown In Figure 5. Although In 
practice, disturbances corresponding to several points In the mission must be 
considered, for this example only the design point corresponding to the maximum 
temperature gradient through the reflector was considered In selecting the actu- 
ator locations. 

The antenna reflector finite element model contains 420 elements which are 
potential sites for actuators. However, It was determined that adequate control 
can be achieved with actuator locations chosen from the 120 sites available on 
the lower surface elements In the reflector. A configuration with an actuator 
at all 120 sites results In value of g ° 0.0016. 

First, we tried to find the best locations for 12 actuators. In an earlier 
study [15], where Intuitive actuator placement was employed, 12 actuators on the 
lower surface gave g = 0.412. The ESPS algorithm yielded g = 0.275, and the 



WOBI algorithm g = 0.3?9, an improvement of 33 percent and 20 percent 
respectively. The scatter of the ESPS method due to choice of initial configu- 
ration was less than 5 percent while the WOBI method produced about 70 percent 
scatter. Convergence typically required fewer than 10 iterations so that fewer 
than 13,000 of the 1.05 x 10^® possible configurations were checked. The 
DeLorenzo algorithm [16] when applied to this problem did not perform well, pro- 
ducing g = 0.459, which is worse than the configuration chosen intuitively in 
[15]. This is likely due to the generally poor performance of this algorithm 
when the number of actuators is much less than the number of possible sites. To 
check whether the result from the DeLorenzo method corresponded to a local mini- 
mum it was used as an initial configuration for the ESPS algorithm. The ESPS 
algorithm moved from the DeLorenzo design and converged to one having g = 0.282 
which is quite close to the best design obtained previously by ESPS. 

The actuator locations from the ESPS, WOBI, and DeLorenzo algorithms are 
shown in Figure 6 along with the locations used in [15]. Figure 7 depicts the 
uncorrected and corrected shapes for the antenna based on the various actuator 
locations. Shown are cross-sections of the shapes corresponding to a section 
through a diagonal of the reflector (the line y = 0). These shapes are intended 
to convey the nature of the shape correction associated with the various actua- 
tor placement algorithms. However, the figure does not include all of the 
points in the structure, and therefore it is not as suitable as the values of g 
for comparing the overall effectiveness of the various algorithms. 

A study was performed to test the behavior of the ESPS, WOBI, and DeLorenzo 
algorithms as the number of actuator locations was increased. The number of 
locations was varied from 12 to 12D out of the 120 available sites. To minimize 
the effect of scatter, the sequential strategy described in the previous section 
was used with the ESPS and WORI algorithms. This strategy was not needed for 
the DeLorenzo algorithm. 
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The results of the comparison of the three algorithms are summarized in 
Table 3. The table shows that the ESPS and WOBI algorithms produced designs 
superior to those obtained by the DeLorenzo algorithm, for all numbers of 
actuators. However, as the number of actuators is increased the DeLorenzo algo- 
rithm requires analyzing a smaller number of configurations than ESPS and WOBI. 
Consequently, when the number of actuators is a large fraction of the number of 
available sites, the DeLorenzo technique would be a reasonably good choice. 

Table 3 suggests that for antenna structures it may be difficult to satisfy high 
surface accuracy requirements with a small number of actuators even if these are 
optimally placed. For example, 60 sites are needed to reduce the distortion by 
an order of magnitude (g * .112 for nc = 60). 

Finally to give a graphical indication of the effect of increasing the num- 
ber of actuators, the corrected shapes for 12 and 40 alctuators located by the 
ESPS technique are shown in Figure B. Use of 40 actuators not only results in 
an increased reduction in overall distortion (g = 0.157 compared to 0.275) but 
also yields a relatively smooth shape compared to the highly oscillatory shape 
produced by 12 actuators. 


CONCLUDING REMARKS 

Two heuristic algorithms were described for the optimal selection of actua- 
tor locations to correct surface distortion of orbiting spacecraft. These algo- 
rithms are denoted Worst-Out-Best-In (WOBI) and Exhaustive-Single Point 
Substitution (ESPS). The algorithms produce results which depend somewhat on 
the intial guess - however, they determine improved locations while evaluating 
only a small fraction of the possible choices. The computational efficiency of 
the algorithms was enhanced by the derivation of fast re-analysis techniques for 
estimating the effect of changing the location of an actuator. 
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The algorithms were demonstrated for a free-free beam and a space antenna 
reflector and the performance of the algorithms was compared to those previously 
obtained with a set of intuitively-located actuators. It was shown that the 
WOBI and ESPS algorithms were able to significantly improve shape corrections by 
relocating actuators. For example, in the beam with 8 actuators the WOBI and 
ESPS corrections were up to 38 percent better than corrections obtained by 
equally-spaced actuators. For the reflector with 12 actuators, WOBI was 20 per- 
cent better than the intuitively-placed actuators and ESPS was 33 percent 
better. As part of the present work, a previously-developed location selection 
technique due to DeLorenzo was also evaluated for the antenna example. While 
the DeLorenzo technique is often computationally cheaper than the WOBI and ESPS 
algorithms, it did not perform as well as either when a relatively small number 
of the available sites had to be selected. When a large fraction of the avail- 
able sites were used, the DeLorenzo technique was a reasonable choice. 
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TABLE 1.- INITIAL DESIGNS USED TO ASSESS EFFECT OF STARTING POINT ON FINAL 
DESIGN FOR BEAM (n = NUMBER OF ACTUATORS, m = NUMBER OF AVAILABLE 
SITES). 


20 


40 


80 


1,2,3.18,19.20 

2,4,6,15,17,19 

3.6.9.12.15.18 

2.3.4.17.18.19 

3.5.7.14.16.18 

3.4.5.16.17.18 
4,5,6,15,16,17 
5,6,7,14,15,16 
6,7,8,13,14,15 
7,8,9,12,13,14 


1,2,3,4,17,18,19,20 

2.4.6.8.13.15.17.19 

2.3.4.5.16.17.18.19 
3,4,5,6,15,16,17,18 

4.5.6.7.14.15.16.17 
8 5,6,7,8,13,14,15,16 

6,7,8,9,12,13,14,15 
7,8,9,10.11 ,12,13,14 

3.5.7.9.12.14.16.18 
4,5,8,10,11,13,15,17 


1,2,3,38,39,40 

2,4,6,35,37,39 

3.6.9.32.35.38 

4.8.12.29.33.37 

5.10.15.26.31.36 
6,12,18,23,29,35 

2.3.4.37.38.39 

3,4,5,36,37,38 

4.5.6.35.36.37 

5.6.7.34.35.36 


1 ,2,3,78,79,80 

2,4,6,75,77.79 

3,6,9,72,75.78 

4,8,12,69,73,77 

5,10,15,66,71,76 

6,12,18,63,69,75 

7,14,21,60,67,74 

8,16,24,54,65,73 

9,18,27,54,63,72 

10,20,30,51,61,71 


1,2,3,4,37.38,39.40 

2.4.6.8.33.35.37.39 

3.6.9.12.29.32.35.38 

4.8.12.16.25.29.33.37 
5,10,15,20,21 ,26,31 ,56 

2.3.4.5.36.37.38.39 

3.4.5.6.35.36.37.38 

4.5.6.7.34.35.36.37 

5,6,7,8,33,34,35,36 
6,7,8,9,32,33,34,35 


1 ,2,3,4,77,78,79,80 
2,4,6,8,73,75,77,79 

3.6.9.12.69.72.75.73 
4,8,12,16,65,69,73,77 
5,10,15,20,61 ,66,71 ,76 
6,12,18,24,57,63,69,75 
7,14,21,28,53,60,67,74 

8.16.24.32.49.57.65.73 
9,18,27,36,45,53,61,69 
10,20,30,40,41 ,51 ,61 ,71 


(SEE FIGURE 1 FOR SITE LOCATIONS) 
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TABLE 2.- SCATTER IN RMS REDUCTION FACTORS g DUE TO STARTING POINT FOR VARIOUS 
NUMBERS OF SITES (m) AND ACTUATORS (n) FOR BEAM PROBLEM. 



n 

m 

ALGORITHM 




WOBI 

ESPS 



20 

0.0508 - 0.0709 

0.0508 - 0.0709 

0.0711 

6 

40 

0.0542 - 0.0843 

0.0499 - 0.0550 



80 

0.0497 - 0.0892 

0.0476 - 0.0544 



20 

0.0358 - 0.0528 

0.0358 - 0.0419 

0.0542 

8 

40 

0.0364 - 0.0601 

0.0341 - 0.0443 



80 

0.0350 - 0.0551 

0.0332 - 0.0383 


TABLE 3.- COMPARISON OF RMS REDUCTION FACTORS (g) AND NUMBER OF CONFIGURATIONS 
EVALUATED (nc) BY THREE HEURISTIC OPTIMIZATION ALGORITHMS FOR 
ANTENNA REFLECTOR. 


NUMBER 

OF ACTUATORS 

12 

15 

20 

30 

40 

60 

80 

120 

ESPS 


0.228 

(42208) 

0.198 

(67978) 

0.179 

(129054) 

0.157 

(248449) 

0.112 

(561876) 

B 

0.0016 

(1) 

WOBI 

0.329 

(1687) 

IB 


0.192 

(6266) 


0.134 

(22047) 


0.0016 

(1) 

DeLorenzo 
(Ref. 16) 

m 



0.249 

(6795) 

0.225 

(6440) 

m 

El 

0.0016 

(1) 

(Ref. 15) 

0.412 









upper entry - g 
lower entry - n^ 

































(a) Beam geometry and nomenclature. 
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12... 19 20 

TWENTY AVAILABLE SITES (m = 20) 

l-/^’2'^3. . . 39^^40 

FORTY AVAILABLE SITES (m = 40) 


zl 


79^^ 80 


. 4 ^ 


EIGHTY AVAILABLE SITES (m = 80) 


(b) Available sites for 20, 40, and 80 actuators. 

Figure 1.- Geometry and locations of available actuator sites for 
beam problem. 
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• * ° ^ (g = 0.0332) 
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I ' (g = 0.0551) 

EQUl SPACED 
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(a) 8 locations chosen 
from 80 sites. 


I o • • 


, ESPS BEST 

(g = 0.0476) 

WOBI WORST 

(g = 0.0843) 


EQUl SPACED 

I (g = 0.0711) 

(b) 6 locations chosen 
from 80 sites. 


Figure 2.- Worst and best actuator locations obtained by WOBI and 
ESPS for beam. 









• DEFORMED SHAPE 

o CORRECTED SHAPE WITH EIGHT EQUALLY- 
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Figure 3.- Deformed and corrected shapes for beam. 



Figure 4.- Tetrahedral truss antenna reflector. 
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Figure 5.- Temperature history for antenna reflector. 


(a) Locations used in 

in reference 15, G = 0.412. 



(b) Locations from ESPS 
method, G = 0.275. 


(c) Locations from WOBI 
method, G = 0.329. 



(d) Locations from DeLorenzo 
algorithm, G = 0.459. 


Figure 6.- Actuator locations for control of antenna surface distortion. 
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Figure 7.- Control of antenna surface distortion by applied temperatures 
(12 actuators). 
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Figure 8.- Corrected shapes with 12 and 40 actuators (locations determined 
by ESf^S). 
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